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1 . INTRODUCTION 


This paper documents a method used for generating uniform chromaticity scale 
imagery from Landsat data. A previous study (ref. 1) was made to map multi- 
channel Landsat data into color space using a maximal chromatic expansion 
method. This study extends the work of reference 1 to the utilization of a 
uniform chromaticity scale (UCS) in the form of color film products. 

Familiarity on the part of the reader with standard colorimetric nomenclature 
is assumed; references 2 and 3 are recommended for the novice. 

The motivation behind generating UCS imagery from Landsat data is uniform, 
controlled perceptibility of color difference caused by differences in 
Landsat data vectors. One can move a certain distance around any color 
center before noticing the difference between the color center and the 
translated point. The distance one can move depends on the direction of 
the motion and on the location of the color center itself. The generally 
ellipsoidal surface surrounding the color center is one step in percepti- 
bility from the color center. The nonlinear transformation to a UCS space 
is such that the ellipsoids become spheres with the same radius at all color 
centers. In this circumstance, the color difference is simply the length of 
the straight line connecting the points whose color difference is desired. 

One of these transformations was used to generate the UCS imagery from the 
Landsat data reported herein. A description of this transformation is given 
in section 2. Formulation of the algorithm used for generating UCS imagery 
is presented in section 3. Conclusions and recommendations of this study 
can be found in section 4. 
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2. DESCRIPTION OF THE UCS TRANSFORMATION USED IN THIS STUDY 

The (L*, a*, b*} UCS transformation (ref. 4) was used in this report. In a 
UCS space, a Euclidean metric is proportional (approximately) to perceptibility 
of color difference. (See reference 2 for a description of some color- 
difference formulas and a more thorough exposition of the subject.) The 
(L*, a*, b*) transformation was adopted by the International Illumination 
Committee [ ( IIC) ; better known as the Commission Internationale de I'Eclairage 
(CIE) ] . 


Let X, Y, and Z be the tri stimulus values; the UCS space is generated from 
the tri stimulus values as follows: 



where Xq, Yq, and Zq define the color of the nominally white object-color 
stimulus; i.e., the illumination spectrum, taken here as that of the light 
table illuminating the transparencies. The L* coordinate may be thought of 
as the lightness or brightness of a color, a* as the green/magenta balance, 
and b* as the yellow/blue balance. 


Slices of the (L* , a*, b*) space (at constant values of a*) that are acces- 
sible to the color gun cube of the production film converter (PFC) were gen- 
erated. Plots of some of those slices are shown in figure 1 for a* values 
of -10, -5, 0, 5, and 10. In these plots, an edge marked B = 0 indicates 
that the blue gun count of the PFC is zero at that edge. An edge marked 
R = 1 indicates that the red gun count of the PFC is 255 at that edge. A 
point in the enclosed area of any of the plots in figure 1 corresponds to 
red (R), blue (B), and green (6) PFC gun counts between 0 and 255. A PFC 

£&JSD1NG page blank not filmed 
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a* = 


10.00 



Figure 1.— Slices of the (L*, 


a*, b*) space at constant values of a* that are accessible to 
the color gun cube of the PFC. 
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Figure 1 Continued 
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Figure 1 Concluded 





product of the (L*, a*, b*) space, sliced at a* = 10, is shown in figure 2. 
Formulas used for calculating gun counts in figure 2 can be found in 
section 3. 

Other UCS transformations exist in the literature {refs. 5, 6). These trans- 
formations might be considered in later research. The (L*, a*, b*) space was 
chosen over those discussed in references 5 and 6 for the following reasons. 

a. The (L*, a*, b*) space has been adopted by the CIE as the principal 
approximation to a true UCS for the object-color solid. 

b. Inversion from the (L*, a*, b*) UCS approximation into the CIE (x, y, Y) 
system is mathematically tractable. According to MacAdam (ref. 7), the 

(£, Jj g) system (ref. 5) is philosophically preferable to the (L*, a*, b*} 
system for this application. However, there are some considerations, such 
as the incorporation of Semmelroth's (ref. 8) crispening factor, that 
create difficulties in using the (£, j, g) system. 
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Figure 2.- A PFC product of the (L*, a*, b*) space sliced at a* = 10. 
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3. FORMULATION OF THE UCS ALGORITHM 

In generating UCS imagery, the light table and film parameters are taken into 
account. The transformation from PFC counts to transmission (wavelength 
dependent) is regarded as stable. The primary colors considered in this 
study are red (R), green (G) , and blue (B). This section is divided into 
three parts. The first part deals with calculating the primaries of light 
table and film chromati cities. The second part deals with fitting a proto- 
type Landsat data structure into the (L*, a*, b*) UCS space. The third part 
presents the algorithm that generates color gun counts from (L*, a*, b*) 
values for PFC product along with four UCS images of a Large Area Crop 
Inventory Experiment (LACIE) segment. 

3.1 COLORIMETRIC DESCRIPTION OF THE TRANSPARENCY GENERATION PROCESS 

The PFC images a cathode-ray tube (CRT) through color filtration onto color 
reversal film, which is then developed. The CRT display is controlled by 
numerical input; for each color, the film density is very nearly linear with 
respect to input counts. The counts-to-density relationship is carefully 
maintained in exposing and developing the film, and the relationship is 
herein regarded as stable. 

An approximation is made that the manner of addition of sequential images 
can be expressed in the form 

X = 2 x. (4) 

i 

where X is the first tristimulus value and i is an index for the red, green, 
and blue PFC primaries. This also exists for Y and Z assuming ,that the 
film system produces the same final result as would be obtained by sepa- 
ration images projectively added. The approach avoids certain real-world 
problems such as masking (ref. 9). (The matter is further discussed in 
reference 10.) More elegant mathematics are available; therefore, a 
follow-on study to investigate the effect of more highly precise work is 
proposed. 
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Measurements of the spectral radiance of the light table and of the wavelength- 
dependent transmission of film patches from full activation of each primary in 
the PFC were made. These measurements were reduced to the CIE chromaticity 
coordinates of the four sources. 


Parameter 

X. 

1 

Light table 

0.3684 

0.4131 

Red 

.5529 

.4458 

Green 

.3603 

.5730 

Blue 

.1919 

.1540 


Under the usual definition (ref. 4) for the reference white, Y, for the light 
table, was taken to be 100. By requiring that the primaries be added at full 
activation to give the light table's own illumination, one finds that the 
tristimulus values for the primaries (fully activated) and for the light 
table are 





/ 89.18 

39.74 

35.90 

13.53 

= I 100.00 

32.04 

57.09 

10.86 

\ 52.89 

.09343 

6.646 

46.13 


(5) 


The activation a. of the i th primary is defined as the ratio of the transmis- 
sion (measured straddling the peak wavelength) at a given input to the trans- 
mission produced by full activation. The tri stimulus values produced by the 
vector activation a are then 



with the 3-by-3 matrix elements taken from equation (5). Conversely, to pro 
duce a given tristimulus value vector, the activation vector 
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3.2 FITTING LANDSAT DATA SPACE INTO THE UCS SPACE 

Two multi temporal LACIE segments were used to generate scatter plots. The 
Landsat data in each uni temporal segment were rotated using the Kauth transfor- 
mation (ref. 11) with a bias. The Kauth brightness K^, Kauth greenness Kg, 
and Kauth yellowness K biases were taken to be 0, 30, and 53 counts, 
respectively. The scatter plots of Landsat data in the (K^, K g ) space for 
segments 1618 and 1645 are shown in figure 3. Figure 4 shows statter plots 
of the data in the (k^j K ) space. The scatter plots in figure 3 were used 
to produce an overall Landsat (K b , K g ) data space scatter envelope. This 
data space scatter envelope contains the raw K b and Kg counts; i.e., the 
bias was removed and was fitted into the (L*, b*) UCS space as shown in 
figure 5. The UCS space, sliced at a* = 10 (figure 2), was used because it 
had the best shape and was near-central ly located. 
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Figure 3.- Scatter plots of Landsat data in the (K b , K ) space for LAC IE segments 1618 and 1645 
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Figure 3.- Concluded. 
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Figure 4.- Continued 
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Figure 4.— Concluded. 




Figure 5.— Fitting Landsat (K^, K^) data space scatter envelope 
into the (L*, b*) UCS space sliced at a* = 10. 


3-18 


The Landsat K y is taken parallel to the a* dimension of the UCS space. Using 
figure 5 and analyzing K g and into two components along b* and L* yields: 

L* * S(Kg sin e + cos e) + A 3 (9) 

b* * S(K g cos e - K b sin 0 ) + A g (10) 

where S is a scale factor, and A 3 and A g are biases. The two points' p 1 and 
p 2 shown in figure 5 were used to calculate e, S, A 3 , and A g in equations (9) 
and (10). Let n = sin 0 ; substituting p^ and p 2 into equations (9) and 
(10) yields: 

(L T - 4> * S "( K gi - K g 2 ) + - \) 01) 

( bf - b £> * S ^(\ - \) * s f b 2 - K b,) < 12 > 


Using equations (11) and (12), 0 is given by: 


0 = tan 


( 


4 - 41( K g, 


-i/ n '. t3n -i \ 1 ‘H g i h) " ( bf ' % 

{/rrj) [fr - ^ + (4 - 4) 


b 2V K b, 





The scale factor S and 

the biases 

A 3 

and A g are 

given by: 

c = . 
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sin e/K 

- k 

, \ + cos 0/ 

k K - K, \ 


l 9 i 

( 

’ 2 ) ( 

b l b 2 ) 

A 3 

■ 4 - s ( 

S 

sin 0 + K. 

b i 

cos ej 

A 6 

- b| + 

K. 

b l 

sin e - K 

4 

cos q \ 

The a* dimension of the 

UCS space 

is 

given by: 



a* = A 7 + S(K y ) 


(14) 

(15) 

(16) 


(17) 
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where A 7 is a bias factor. The overall scatter plot of figure 4 shows that 
the data center is at raw K y of -8 counts. This center should correspond to 
a* of 10 where the UCS space was sliced to fit the Landsat data as shown in 
figure 5. Hence, A 7 is given by: 


Let 


A ? = 10 + 8S 




= S 



cos 0 
-sin 0 


08 } 


09) 


then, fitting Landsat (K^, K^, K^) data space into the UCS space is achieved 
through the following affine transformation: 

u ~ Fe^ + (20) 

where 


U[ = a column vector of L*, b*, and a* 
e = [K K, K ]', and e‘ is the transpose of e_. 

y w j 

— = ^ A 3 A 6 A 7 ^ ' 

The matrix F is defined as: 



( 21 ) 


Using the numerical values of p-j and p£ of figure 5, F and 6. are computed to 
be: 

( -0.3012 1.0267 0 
1.0267 .3012 0 

0 0 1.0670 
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M 1 . 594' 


6 = 1 -32.288 
1 8.560> 


( 23 ) 


3.3 ALGORITHM THAT GENERATES COLOR GUN COUNTS FROM (L*, a*, b*) PICTURE 
ELEMENT (PIXEL) VALUES FOR PFC PRODUCT 

It was shown in section 3.2 that given the e_ vector for a picture element 
(pixel) in a Landsat scene, the L*, b*, and a* values for that pixel were 
computed according to equation (20). The tri stimulus values X, Y, and Z for 
a pixel are computed from the pixel's L*, b*, and a* values by inverting 
equations (1) through (3) as follows: 


y » Y «-* * 16) 3 
1 l 0 1562500 


(24) 



Z = 




(25) 


(26) 


Let a.. , where i = R, G, and B primaries, be activations proportional to the 
tristimulus values. They are proportional to radiance and, hence, to film 
transmission t. Film transmission and optical density D are related by: 


or 



D = - lo 9l0 



(27) 


The activations a^ (0 < a^ < 1) are produced from the tristimulus values of a 
pixel's color according to the following: 

a_ « T _1 t (28) 
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where 


a [a R a fi a B ]‘ 
t = [X Y Z] 1 

( 39.74 35.90 13.53 

32.04 57.09 10.86 

.09343 6.646 46.13, 

The' film transmission and the pixel's activations are linearly related by: 

t. = x. . + a. (t. - t. . ) ; i = R, G s and B (29) 

l i,min i v i,max i,min' * ’ ’ 


where r- and t- were measured and found to be: 
i ,min i ,max 

T R, max = °-l 762202 ; TR , min = 0.0030142 ; t S i1BX - 0.4933794 ; 

T S. nrin = °- 0065697 ; T B, m ax ‘ 0 - 2020853 • T B,n»1n = °- 0030685 


(30) 


Let 


*^i ,max ~^ og 10 T i ,min 
AD i ” log 10^ T i ,max/ T i ,min^. 


i = R, G, and B 


(31) 


then the color gun counts C.. , i = R> G, and B of the PFC for the pixel, are 
given by: 

V 255\ 


C i = ( log 10 T i + D i ,max)( AD i ) ; i = R, G, and B 


(32) 


Equations (20), (24), (25), (26), (28), (29), and (32) constitute an algorithm 
for generating UCS imagery from Landsat data. These equations were used to 
produce the UCS color film products shown in figure 6 for LACIE segment 1618. 
Figure 6 shows four images of segment 1618, corresponding to four multitemporal 
acquisitions and covering all the biological stages of small grains. Figure 3, 
with figure 5, provides an interpretive key for the images in figure 6; i.e., 
small-grains fields in biostage 1 appear to have a blue color [low greenness 
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Figure 6.- UCS Film product of LACIE segment 1618. 
Acquisition date: May 6, 1976. 
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Figure 6.- Continued. Acquisition date: June 17, 1976. 
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Figure 6.- Continued. Acquisition date: July 17, 1976. 
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Figure 6.— Concluded. 


Acquisition date: Aug. 22, 1976. 
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value Kg j while they change from light yellow orange |Kg^ > K g ^j to a darker 

yellow orange / K > K \ in going from biostage 2 to 3. In biostage 4 
\ ^3 ^ 2 / 

IK < K the small-grains fields appear to have a pinkish color. This is 
\ 9/1 9 -? 


l 94 93 J 

in agreement with how the data were fit into the UCS space as shown in figure 5 


and given by equations (9), (10), and (17). The pinkish cast shows an increase 
in the Kauth yellow component; K is scaled toward positive a*, which is in the 

J 

general direction of the PFC red gun. 
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4. CONCLUSIONS AND RECOMMENDATIONS 


An algorithm for generating uniform chromaticity scale imagery from Landsat 
data has been presented. A computer program was written to implement the 
algorithm, and UCS film products were generated. The colors in the film and 
their temporal change are consistent with what was expected for the particular 
scaling of Kauth components into the (L*, a*, b*) color space. The UCS film 
product has yet to be put to the practical test of competing with previous 
transformations. In that competition {to be done outside the purview of 
this report), the philosophically satisfying notion of transforming Landsat 
data so that a one-count difference is equally perceptible at all locations 
in data space will be tested. 

The authors recommend that analyst-interpreters (AI's) test the UCS imagery 
using a variety of LACIE segments. Preliminary examination indicates that 
the UCS product offers the possibilities of the following. 

a. A single film product that will supplant two film products in current use. 

b. Improved visibility of data differences in regions in data space which are 
critical to crop identification. 

c. An analytic route to the determination of data-space transformations that 
will be optimal for particular discrimination problems. For example, in 
another project, the transformation has been used to display water bodies 
in Landsat data, with encouraging results. 
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